Probing the BCS-BEC crossover with photons in a nonUnear optical fiber 
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We propose a scheme where strongly correlated photons generated inside a hollow-core one- 
dimensional fiber filled with two atomic species can be used to simulate the BCS-BEC crossover. 
We first show how stationary light-matter excitations (polaritons) in the system can realize an 
optically tunable two component Bose-Hubbard model, and then analyze the optical parameters 
regime necessary to generate an effective Fermi-Hubbard model of photons exhibiting Cooper pair- 
ing. The characteristic correlated phases of the system can be efficiently observed due to the in situ 
accessibility of the photon correlations with standard optical technology. 



Superconductivity is undoubtedly one of the most fas- 
cinating and elusive condensed matter phenomena yp. 
Bardeen- Cooper- Schrieffer (BCS) theory [2] provided the 
first satisfying explanation of the effect, by proposing 
that fermions form long range pairs (Cooper pairs) un- 
der an arbitrarily weak attractive interaction. A min- 
imal model exhibiting Cooper pairing is the attractive 
Fermi-Hubbard (FH) model [3, 4 . In the BCS-like re- 
gion characterized by weak inter-species attraction, large 
cooper-pair-like states form the ground state and exhibit 
long-range correlations. The latter is destroyed as the 
inter-species attraction is raised and localized bosonic 
molecules are formed. Seminal experimental realizations 
of FH model require dilute fermi gases for temperatures 
well below the degeneration temperature, making these 
experiments extremely challenging [5 . 

A different approach involves utilizing the well known 
mapping of one-dimensional hard-core bosons into free 
spinless fermions [6 . This led to works suggesting 
bosonic mixtures on an optical lattice in the regime of 
strong intra-species repulsion for an effective realization 
of the necessary interacting fermionic behaviour [7]. In 
this case, the so-called BCS-BEC crossover could be ob- 
served for a higher temperature and furthermore a new 
phase appears as the inter-species attraction is increased 
in comparison to the intra-species repulsion. In the latter 
regime, the system moves away from the fermionic BCS- 
BEC regime and enters a new strongly localized bosonic 
phase termed big boson (BB) with almost all the bosonic 
molecules occupying the same site. 

In this work we show that one could circumvent the is- 
sues around bosonic or fermionic atoms and actually use 
photons in a nonlinear optical medium to efficiently sim- 
ulate the crossover. We show how to generate a highly 
tunable two component Bose-Hubbard (BH) model of po- 
laritons or photons in a nonlinear fiber and analyze the 
probing of BCS-BEC crossover using optical methods. 
Our work is motivated by recent progress in the field of 
photonic quantum simulations [8 , where photon crystal- 
lization, Luttinger liquids behaviour and the "Pinning 
transition" has been recently predicted in hollow core 
optical fibers [9]. 
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Figure 1: (a) A schematic diagram of the system under study. 
A hollow-core optical waveguide is filled with a mixture of 
two species of cold atoms interacting with a pair of quantum, 
^(1,2), and a pair of classical, r^(a,6), light fields. The resulting 
stationary light-matter excitations in the waveguide can be 
steered to a strongly interacting regime mimicking an effective 
Fermi-Hubbard model with highly tunable attractive interac- 
tions, (b) The atomic level structure, (c) Coherently mapping 
the stationary excitations to propagating photon pulses allows 
for the efficient probing of the BCS-BEC crossover by mea- 
suring the temporal correlations of the photon pulses leaving 
the fiber. 

In the following, we introduce a possible experimental 
setup and analyze the conditions for the realization of a 
two component BH model of photons. We then investi- 
gate the possibility of tuning the photonic intra-species 
repulsion to the necessary regime to generate effective 
fermionic behaviour with attractive inter-species inter- 
actions. Our photonic system is then shown to exhibit 
the BB-BEC-BCS crossover as the optical parameters 
are varied within realistic regimes. Finally, we discuss 
a possible observation of the crossover by analyzing spe- 
cific correlation functions of the released photon pulses 
as they exit the fiber. 

The system: As shown in Fig. [l] we consider a waveg- 
uide filled with two species, a and 6, of cold atoms; e.g.^ 
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a hollow-core photonic crystal fiber doped with Rubid- 
ium isotopes [TO^. Our scheme involves four hyperfine 
levels for each atomic specie, two pulsed quantum fields, 
^(1,2) (^7^)7 and two classical laser fields, r^(a,6) (t). 

Initially, the medium is illuminated by the leftward 
propagating quantum pulses £^(1,2),+ (^7 ^) and classical 
fields l](ci,6),+ (^); henceforth, the plus (minus) subscript 
will be used to represent leftward (rightward) propa- 
gation. After the quantum pulses £^(1,2),+ (^7 ^) com- 
pletely enter the medium, switching off the classical fields 
^(a,6),+ (t) adiabatically converts the quantum pulses 
into atomic excitations in the usual slow light manner 
[TTJ[T2]. Then, simultaneously switching on the four clas- 
sical fields ^(a,b),± (t) creates a Bragg grating that traps 
the quantum fields [11]. In the following, we show how 
the trapped fields experience strong Kerr nonlinearities 
induced via their interaction with the atoms, steering 
the system to a strongly interacting regime described by 
a two component Lieb-Liniger (LL) model. 

In the interaction picture, the Hamiltonian of our sys- 
tem shown in Fig. 1 is 
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with 
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where x = a^b denote different atomic species. The con- 
tinuous collective atomic spin operations, = a^^ (z, t), 
give the average of \p)^ {q\ over the x-type atoms in a 
small but macroscopic region around spatial coordinate 
z. The wavevectors are given by k^u'^\ k^^'^\ and their 



central frequencies by, cjqu' , uJ^i' , for the slowly vary- 
ing quantum fields and classical fields, in that order. For 
notational simplicity, we omit the space and time de- 
pendence of the operators. We denote the atomic den- 
sities and denotes the single-photon-single-atom 
coupling constant between an x-type atom and the jth 
quantum field. Here, we have assumed that the quan- 
tum fields ^(1,2) drive the transitions |2)^ (1| and |4)^ (3| 
with the same strength g^-^ 2) • The one-photon detunings 
are denoted as Af and Af; the two-photon detunings as 
Af; the quantum pulse detunings are written as Sj with 



= uj^J - uj^^J, Si = -S^, and (^f = , 
Similarly to earlier slow-light schemes [12], we de- 
fine polariton operators ^j,± = g^^ V27rn^Ej^±/Qx.^ 
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where Xi^2 = cl-,^-, that describe the long lived sta- 
tionary light-matter excitations of the system. Insert- 
ing these into the relevant Maxwell-Bloch equations de- 
scribing the propagation of the quantum fields in the 
medium, we find that the stationary polariton opera- 
tors, ^1^2 = (^(1,2),+ + ^(i,2),-)/2, obey two coupled 
nonlinear Schrodinger equations originating from a two- 
component Lieb-Liniger Hamiltonian: 
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V}^^ cos^ (nni^^ z^ + Xj ^] ^] } 

-FXi2 J (i2;^l^i^5^2. 

Here, the effective masses of polaritons are mj = 
-Au^^^ /{2u^^^vl^^)-T%n^^^^ /{AA^'vl^^) with Acj^^'^) = 

u^/^ - c^r'^^ and 4^'^^ = -''^''^U,^^^^^^^ 
are the group velocities of the quantum fields in the 
medium with the corresponding velocities in an 

empty waveguide. The spontaneous emission rate 
of a single atom into the waveguide modes ^ii) — 
A7i{g^-^Y / . The intra-species repulsions are Xj = 

(A^^ )^S^j r^^'U^"^^/(2A4^ ) and the interspecies repulsions 



IS X12 - 2n^{/s.l-5^magl) 
^Xj _ o2 l(n'2 
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^V^Vm and S-.- = (A^^- - 
A3^/2)/(A4^ — Ag^). The Hamiltonian above contains 
a periodic potential term induced by a slight modulation 
in the atomic density distributions for the two atomic 
species [9 . This is possible, for example, by applying 
two external fields. Thus the atoms, initially prepared 
in their ground states |1)^^'^\ have densities given by 



cos'^(7rni^^^z), where 



Ax 3) 



is the linear 



density of sites for atomic specie Xj, and n^^ ^ n^^ ; 
ie., the modulation is only a perturbation. Such peri- 
odic distributions for the atomic species give the trapping 

polaritonic potential Vq 



A3) 
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and the periodic polaritonic lattice with depth V{ 
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A two component Bose-Hubbard model of polaritons: 
For a sufficiently strong periodic potential and weak in- 
teraction between the polaritons, the LL Hamiltonian 
above can be mapped to a two species Bose-Hubbard 
model. For simplicity, we assume Acj^-^^ = and that 
the counter propagating classical fields are identical, and 
the two atomic species have identical mass, distribution, 
and interaction with the corresponding quantum fields, 
U) _ ^, o — o oo^-^') 



^1' ^IB = TiD, 



g and A^ 



Ak for 



k = 2,3,4. The conditions on the optical parameters 
regime for the mapping from LL to BH to be valid trans- 
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late to the single photon detuning /S^/^/V > 20 and control 
laser Rabi frequencies of l^/F < 3. T is the atomic decay 
rate which is assumed at 20MHz for the typical Rb tran- 
sition in question. The above conditions are calculated 
following the methods in [9], for atomic densities equal 
to n = 10^m~^ with 10% modulation ni = O-ln^; pho- 
tonic densities of = 10^ m~^, single atom cooperativity 
T] = Fid/F = 0.2 and two photon detunings A2 = — 5F, 
A3 = -O.OIF [13] 

Assuming that the two-photon detuning obeys IA3I <C 

^ , A4, SO that the dimensionless quantities A-^ , S-^ 1, 
we obtain a two component BH model for polaritons: 



H 
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■vY^it^a^ (4) 



where a^^r is an annihilation operator of a a-type po- 
lariton at ith site and stands for nearest neigh- 

bors, t = W^^^E]l^exp (^-2y^Vi/ Er^ /v^ = = 

t;, U = V2^xins{Vi/ER)^^^ = = U^, V = 
V27rxi2ns{Vi/ER)^^^ /2, where Er = 7r^nl/{2m) with 
m = — FiDn/ {4:A2Vg) and = density of sites; x = 



FiD^^A. 



TiBVg/ (2A4) = xi = X2, X12 ^ 
andFi = -FiDA3^^ni/(41]2). 

To discuss the simulation of BCS-BEC crossover, we 
focus on the case of repulsive intra-species {U > 0), 
and attractive inter-species (V < 0) interactions, which 
can be achieved by setting x > and Xi2 < 0- The 
ratios between inter- and intra-species interaction V/U 

= X12/X1 = a2 Aa^2 and the hopping to repulsion ratio 

t/U = 4Fi'/'^^/'exp (^-2y/W/E^) I (V2^Xin.) deter- 
mine the physics of the Hamiltonian Q completely. As- 
suming the effective mass to be positive, the constraints 
on the one-photon detunings to probe the BCS regime 
are A2 < and < A4 < (^2- We note here that dif- 
ferent regimes are also possible that can lead to different 
effects such as spin-charge separation or Kondo physics. 

From bosonic polaritons to paired fermions: Figure 
[2] shows the accessible regimes of the inter- and intra- 
species interactions such that the BH model of polaritons 
can be mapped into an effective FH model. It is shown 
that by simply tuning quantum optical parameters such 
as the single photon detunings A4, the detuning between 
the quantum fields (^2 ? strength of the classical trap- 
ping lasers 1^, the required regime, t/U <C 1, can be 
reached, while leaving a range of values of V/U accessi- 
ble. 

Under these restrictions, the polaritonic two-species 
Bose-Hubbard model shows the fermion-like BEC-BCS 
crossover as well as the highly bosonic BB behaviour 
[7]. Preparing these BCS- and BEC-like states is possible 
here by simple tuning of the optical parameters, forcing 
the trapped polaritons to form large Cooper-pair-like ob- 
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Figure 2: (a) Inter-species interaction ratio V/U and (b) hop- 
ping parameter ratio t/U, characterizing the realizable two 
component BH model of photons as a function of the quan- 
tum pulse S2 and single photon A4 detunings. Fid = 0.2F, 



i/us = 10^, n/m = 10, A2 = -5F, and A3 
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jects or localized bosonic molecules, on demand. Switch- 
ing off the trapping lasers and coherently mapping the 
stationary polaritonic correlations to propagating pho- 
ton pulses allows one to probe the different states by 
observing the temporal (and hence spatial) second-order 
correlation functions of the leaving the waveguide. We 
note here the in situ character of the measuring process 
in contrast to the cold atom cases where, usually, time of 
flight imaging is needed [3 . 

To discuss the crossover, we focus on the second or- 
der cross-species correlations (l) = and 
the correlations between cross-species population differ- 

(2) 

ences g_^{l) = Ei((^it " - ^i+n)) as func- 

tions of the site number difference The dependence 
on t/U and V/U is implied. These types of correlation 
functions can be measured by collecting the component- 
resolved photon-counting records and analyzing the col- 
lected data. For example, one could use a beam splitter 
and energy-resolving photon detectors to collect the re- 
quired data. The strength of our proposal lies in the pos- 
sibility of directly observing the spatial correlations that 
are the defining characteristics of the different phases. 

To calculate the correlation functions, the ground state 
of Eq. (4) is computed numerically for 6 polaritons in 8 
sites. This corresponds to assuming 3 photons in each 
quantum field entering the waveguide and a polaritonic 
potential modulation with the wave vector km = 27r*8/L, 
with L the length of the fiber which is taken to be a 
few centimeters. The BB-BEC-BCS crossover can eas- 

f2) 

ily be seen from the on-site correlation ^|^^(0) as shown 

r2) 

on Fig. 3. Abrupt changes in ^^^^(0) (normalized to the 
value at |V^|//7 = 1.5) at \V\ = U indicate a transition 
from the big boson (BB) state, where all the polaritons 
pair up at a single site, to the localized pairing (BEC) 
state, where different pairs prefer to space out. The 
curves also indicate a crossover from the locally paired 
(BEC) state (when |V| ^ t) to the long-range paired 
(BCS) state (|V| <C t). The coloured background por- 
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Figure 3: Cross-species correlations (/) and population dif- 
ference correlations g^^\l) for (a) t = 0.01 and (b) t = 0.1 as a 
function of the inter-species interaction calculated at the same 
site / = and for neighbouring sites / = 1. The colored gra- 
dient background, proportionate to the on-site cross-species 
correlation, portrays the BB-BEC-BCS crossover. By map- 
ping the polaritons to photons and releasing the pulses, the 
different phases of the system can be probed using photon 
intensity correlation measurements. 
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Figure 4: Cross-species second order correlation function for 
t/U = 0.01 and different values of \V\/U = 1.2,0.99,0.5,0.01 
as a function of the distance in units of the effective photonic 
lattice spacing / . 



Conclusion: We have shown that slow-light-EIT-based 
techniques can be used to generate a two-species ID 
Bose-Hubbard model of photons inside a hollow-core 
waveguide filled with cold atoms. The tunability of 
the optical parameters allows for the strongly repulsive 
regime to be reached, where the predictions of the ID 
Fermi-Hubbard model can be studied experimentally. In 
particular, the resulting strongly correlated polaritons 
(or photons) in the waveguide can mimic the BB-BEC- 
BCS crossover by appropriate tuning of the optical pa- 
rameters. The possibility of in situ measurements of 
the spatial correlation functions makes our scheme an 
excellent candidate for an efficient observation of the 
crossovers using light. 



trays the different phases and how they cross over. The 
black dashed curves show the correlations at / = 1, in- 
dicating that the photonic BEC pairs space out in the 
BEC region as is increased until the system crosses 
over to the BCS regime. The white dotted and dash- 
dotted curves illustrate 5'^_)(0, 1) normalized to the value 
at |V^|//7 0. They show high sensitivity to the BCS- 
BEC crossover, but are not suitable for observing the 
BB-BEC crossover. In the in situ photon correlation 
measurements, the BB-BEC-BCS crossover will appear 
as a transition from a highly bunched behaviour in the 
BB regime, to strongly anti-bunched behaviour in the 
BSC regime in the on-site cross-species intensity correla- 
tions g[l\0). 

While the correlation functions at / = 0, 1 give a good 
signpost for the three phases, correlations at longer dis- 
tances are required, especially in the BCS regime, to com- 
pletely describe the physics. Figure 4 shows (l) for 
different values of \V\/U with t/U = WM fixed. The ex- 
pected short-range to long-range crossovers are clearly 
visible, indicating the measurable signs of the elusive 
BCS-BEC crossover in our scheme. 
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